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ABSTRACT 

We present a theoretical model for supernova (SN) 2008D associated with the luminous X-ray tran- 
sient 080109. The bolometric light curve and optical spectra of the SN are modelled based on the 
progenitor models and the explosion models obtained from hydrodynamic/nucleosynthetic calcula- 
tions. We find that SN 2008D is a more energetic explosion than normal core-collapse supernovae, 
with an ejecta mass of M^j = 5.3 ±1.0 Mq and a kinetic energy of Ek = 6.0 ± 2.5 x 10^^ erg. The pro- 
genitor star of the SN has a 6 — 8Mq He core with essentially no H envelope (< 5 x 10~^ Mq) prior to 
the explosion. The main-sequence mass of the progenitor is estimated to be Mms = 20 — 25 Mq , with 
additional systematic uncertainties due to convection, mass loss, rotation, and binary effects. These 
properties are intermediate between those of normal SNe and hypernovae associated with gamma-ray 
bursts. The mass of the central remnant is estimated as 1.6 — 1.8Af0, which is near the boundary 
between neutron star and black hole formation. 

Subject headings: supernovae: general — supernovae: individual (SN 2008D) — nuclear reactions, 
nucleosynthesis, abundances — radiative transfer 



1. INTRODUCTION 

A luminous X-ray transient was discovered in NGC 
2770 in the Swift XRT data taken on 2008 January 9 for 
the observation of SN 2007uy in the same galaxy (Berger 
& Soderberg 2008; Kong & Maccarone 2008). The X- 
ray emission of the transient reached a peak ^ 65 sec- 
onds, lasting ~ 600 seconds, after the observation started 
(Page et al. 2008). The X-ray spectrum is soft, and no 
7-ray counterpart was detected by the Swift BAT (Page 
et al. 2008). 
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Given the small total X-ray energy and the soft X- 
ray emission, Soderberg et al. (2008) and Chevalier & 
Fransson (2008) interpreted the X-ray transient as a su- 
pernova (SN) shock breakout. On the other hand, Xu et 
al. (2008), Li (2008) and Mazzah et al. (2008) considered 
this transient as the least energetic end of gamma-ray 
bursts (GRBs) and X-ray flashes (XRFs). 

The optical counterpart was discovered at the posi- 
tion of the X-ray transient (Deng & Zhu 2008; Valenti et 
al. 2008a), confirming the presence of a SN, named SN 
2008D (Li & Filippenko 2008). To study this event in 
detail, intensive follow-up observations were carried out 
over a wide wavelength range including X-rays (Soder- 
berg et al. 2008; Modjaz et al. 2008b), optical/NIR 
(Soderberg et al. 2008; Malesani et al. 2009; Modjaz et 
al. 2008b; Mazzah et al. 2008), and radio (Soderberg et 
al. 2008). 

SN 2008D showed a broad-line optical spectrum at 
early epochs (t < 10 days, hereafter t denotes time after 
the transient, 2008 Jan 9.56 UT; Soderberg et al. 2008). 
However, the spectrum changed to that of normal Type 
lb SN, i.e., SN with He absorption lines and without H 
lines (Modjaz et al. 2008c). To date, the SN associated 
with GRBs or XRFs are all Type Ic, i.e., SNe without H 
and He absorption. 

Soderberg et al. (2008) and Mazzah et al. (2008) 
estimated the ejected mass and the kinetic energy of 
SN 2008D. Using analytic formulae, Soderberg et al. 
(2008) suggested that this SN has the ejecta mass Mgj = 
3 — 5Mq and the kinetic energy of the ejecta Ek ~ 
2 - 4 X 10^1 erg. Mazzah et al. (2008) did model cal- 
culations and suggested that this event is intermediate 
between normal SNe and GRB-associated SNe (or hyper- 
novae), with A/oj ~ 7Af0 and Ek ^ 6 x 10^^ erg. In their 
model calculations, hydrodynamic/nucleosynthetic mod- 
els are not used, and thus, their estimate of the core mass 
prior to the explosion and the progenitor main-sequence 
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^Main-sequence mass (Mq) estimated from the approximate formula obtained by Sugimoto & Nomoto (1980) 

•^The mass of the He star (Mq) 

•^Progenitor radius prior to the explosion (Rq) 

•^Mass cut (Mq) 

•^The mass of the SN ejecta (Mq) 

The kinetic energy of the SN ejecta (10^^^ erg) 
SThe mass of ejected ^^m (Mq) 
''Velocity at the bottom of the He layer (km s"^) 
'Velocity at the outer boundary of ^^Ni distribution (km s^^) 
JThe ejecta mass at u > 0.1c (Mq) 
''Estimated from the photospheric radius and temperature of the early part of the LC (i < 4 days) using the formulae by Waxman et al. 

(2007). This is consistent with the estimate by Modjaz et al. (2008b) while Chevalier & Fransson (2008) derived ^^ 9Rq using Af^j and 
i?K estimated by Soderberg et al. (2008) and the formulae by Chevalier (1992). 



mass is less direct. 

In this paper, we present detailed theoretical study 
of emission from SN 2008D. The bolometric light curve 
(LC) and optical spectra are modelled based on the pro- 
genitor models and the explosion models obtained from 
hydrodynamic/nucleosynthetic calculations. In §2, we 
show the progenitor and explosion models. The opti- 
cal LC and spectra are modelled in §3 and §4, respec- 
tively. We discuss the nature of SN 2008D in §5 and 
finally give conclusions in §6. Throughout this paper, we 
adopt 31 Mpc {fi = 32.46 mag) for the distance to SN 
2008D (Modjaz et al. 2008b; Mazzali et al. 2008) and 
E{B — V) ^ 0.65 mag for the total reddening (Mazzali 
et al. 2008). 

2. MODELS 

To understand the nature of SN 2008D and its pro- 
genitor star, (1) we first construct the exploding He 
star models with various masses by performing hydrody- 
namic/nucleosynthetic calculations for the presupernova 
He star models. (2) The important parameters of the SN, 
such as the mass (Mej) and kinetic energy of the ejecta 
(£'k), are estimated by modelling the bolometric LC (§3) 
and the optical spectra (§4). Then, (3) the He star mass 
(Mq) can be estimated from the best set of Mej and Ek- 
Finally, (4) the main-sequence mass (Mms) of the pro- 
genitor star can be estimated by the evolution models, 
which predict the Mus-Ma relation. 

In this section, we construct hydrodynamic models 
using an evolutionary model. The progenitor model 
and hydrodynamic/nucleosynthetic calculations are de- 
scribed in §2.1 and §2.2, respectively. The hydrodynamic 
models are tested in §3 and §4. 

2.1. Progenitor Models 

In the strategy described above, presupernova He star 
models are required as input for the hydrodynamic cal- 
culations. The Wolf-Rayet star models with stellar winds 
tend to form the He stars whose masses are larger than 
that inferred for SN 2008D (e.g., Soderberg et al. 2008). 



To study the properties of the ejecta and the progen- 
itor star without specifying the mass loss mechanism 
(stellar winds in Wolf-Rayet star or Roche lobe overflow 
in close binary, i.e., possible binary progenitor scenario, 
e.g., Wellstein & Langer 1999), we adopt the He star 
evolution models with various masses. We use five He 
star models with the masses of 4, 6, 8, 10 and I6M0 
(Nomoto & Hashimoto 1988; Nomoto et al. 1997; Naka- 
mura ct al. 2001b), where the mass loss is not taken 
into account. These models are called HE4, HE6, HE8, 
HEIO and HE16, respectively. The corresponding main- 
sequence masses of these models are w 15, 20, 25, 30 
and 4OM0, respectively (Table 1), which is estimated 
from the approximate formula of the Mms-M^ relation 
obtained by Sugimoto & Nomoto (1980; Eq. 4.1). 

The difference in the density structure of the He core 
is negligible among different stellar evolutionary calcu- 
lations if the He core mass is the same. As a result, 
the observable quantities (i.e., LC and spectra) after the 
hydrodynamic and radiative transfer calculations are not 
affected by the variety of the evolutionary models. Thus, 
the estimate of Mej and Ek does not depend on the evo- 
lutionary models. 

However, the Mms-^^q relation depends on the several 
evolutionary processes that are subject to uncertainties, 
e.g., the convective overshooting, wind mass loss, shear 
mixing and meridional circulation in rotating stars. The 
different assumptions adopted in different stellar calcu- 
lation codes may affect the Mms-Mqi relation as summa- 
rized in Figure 1 . The red line shows the formula derived 
by Sugimoto & Nomoto (1980). The other lines show 
the relations obtained by the models including mass loss 
(Limongi & Chiefh 2006, blue; Hirschi et al. 2004, green; 
Rauschcr ct al. 2002, black), and rotation {vzams = 300 
km s^^; Hirschi et al. 2004, cyan). The models shown 
in open squares have a H envelope prior to the explosion 
while the models shown in the filled squares have a bare 
He core. 

These models assume the solar abundance for the ini- 
tial abundance. The mass loss causes the smaller M„ 
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Fig. 1. — Relation between the main-sequence mass and the He 
core mass in the formula derived by Sugimoto & Nomoto (1980; 
red, used in Nomoto & Hashimoto 1988). It is compared with the 
relation between the main-sequence mass and the He core mass at 
the presupernova stage in the models by Limongi & Chieffi (2006, 
blue), Hirschi et al. (2004, green and cyan for non-rotating and 
rotating models, respectively), and Rauscher et al. (2002, black). 
The models shown in open squares have H envelope while the mod- 
els shown in filled squares have a bare He core. 
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Fig. 2. — Density profile of the explosion models at one day 
after the explosion. Red (HE4), blue (HE6), green (HE8), magenta 
(HEIO), and gray (HE16) lines show the models for SN 2008D (see 
Table 1). The black lino shows the C+O star explosion model used 
for SN 1998bw in Nakamura et al. (2001a). The vertical lines show 
the velocity at the bottom of the He layer in each model (with the 
same colors). 



for Mms > 30Mq in the models by Limongi & Chieffi 
(2006) and Hhschi et al. (2004). For the lower metal- 
licity models, the mass loss rate is lower, so that the 
Mus-Ma relation would be closer to that of Sugimoto 
& Nomoto (1980). For the stars with A/ms < SOAfg or 
Ma < IOMq, the gradient in the plot is similar among 
the models. 

At the presupernova stage, the He stars consist of the 
Fe core, Si-rich layer, 0-rich layer, and He-rich layer. 
The mass of the Fe-core is ^ 1.4 — 1.6AfQ, depending 
on the model. The mass of the 0-rich layer is sensitive 
to the progenitor mass, while the mass of the He-rich 
layer is ^ 2Mq irrespectively of the He star mass. Note 
that the mass of He-rich layer can be as large as ~ 3Mq 
depending on the evolutionary models (e.g., Limongi & 
Chieffi 2006) and can also be smaller than 2Mq prior to 
the explosion by mass loss. 

The mass fraction of O in the 0-rich layer is ^ 0.8. 
Other abundant elements in this layer are Ne, Mg, and 
C, with mass fractions of order 0.1. These are almost 
irrespective of the evolutionary models. The He mass 
fraction in the He-rich layer is ~ 0.9. The second most 
abundant element in this layer is C, with a mass fraction 
of ~ 0.03, but this is rather uncertain (§3). Oxygen is 
also produced in the He-rich layer, but the mass fraction 
of O is only -0.01. 

2.2. Hydrodynamics & Nucleosynthesis 

The hydrodynamics of the SN explosion and explo- 
sive nucleosynthesis arc calculated for the five progenitor 
models. The hydrodynamic calculations are performed 
by a spherical Lagrangian hydrodynamic code with the 
piecewise parabolic method (PPM; Colella & Woodward 
1984). The code includes nuclear energy production from 
the a network. The equation of state includes gas, radia- 
tion, e~-e^ pairs. Coulomb interaction between ions and 
electrons and phase transition (Nomoto 1982; Nomoto & 
Hashimoto 1988). The explosion is initiated by increas- 



ing the temperatures at a few meshes below the mass cut 
(see below), i.e., a thermal bomb. 

The SN ejecta become homologous at — 1000 s af- 
ter the explosion. After the hydrodynamic calculations, 
nucleosynthesis is calculated for each model as a post- 
processing (Hix & Thielemann 1996, 1999). The reaction 
network includes 280 isotopes up to ^^Br. The results of 
the nucleosynthesis depends on the progenitor mass and 
the kinetic energy of the explosion. The kinetic energies 
in the five models are determined to explain the observed 
LC (§3). 

The explosion models are summarized in Table 1. The 
mass cut (Afcut) is defined after the nucleosynthesis cal- 
culation to eject the optimal amount of '^^Ni to power the 
LC. Figure 2 shows the density structure of the explosion 
models at one day after the initiation of the explosion. 
The "bump" in the density profile is caused by the re- 
verse shock generated at the boundary of the C+O/He 
layers. 

The vertical lines in Figure 2 show the velocity at the 
bottom of the He layer after the expansion of SN ejecta 
become homologous (wHe, Table 1). Since strong mixing 
is expected in less massive stars (Mq, < AMq, Hachisu 
et al. 1991), the value of ^He in model HE4 is the up- 
per limit of the inner velocity of the He-rich layer. If 
the mass of the He layer prior to the explosion is larger 
(smaller) than ~ 2Mq as in our model set, vue can be 
lower (higher). 

3. BOLOMETRIC LIGHT CURVE 

The pseudo-bolomctric (UBVRIJHK) LC was con- 
structed by Minczaki et al. (in preparation, see Ap- 
pendix A) compiling optical data taken by the MAG- 
NUM telescope (Yoshii 2002; Yoshii, Kobayashi & 
Minczaki 2003), the Himalayan Chandra Telescope, and 
Swift UVOT ((/-band, Soderberg et al. 2008), and also 
NIR data taken by the MAGNUM telescope. The first 
part of the LC (i < 4 days) seems to be related to the 
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Fig. 3. — Pseudo-bolomctric (UBVRIJHK) LC of SN 2008D (Minczaki ct al. in preparation) compared with the results of LC 
calculations with the models HE4 (red), HE6(blue), HE8 (green), HEIO (magenta) and HE16 (gray). The pseudo-bolometric LC is shown 
in filled (left) and open (right) circles. The thin black line shows the decay energy from ^®Ni and ^®Co [A'/(^^Ni) = 0.07 Mq\. The 
bolometric magnitude at f ~ 4 days after the X-ray transient is brighter by ~ 0.25 mag than that shown by other papers (Soderberg ct al. 
2008; Malesani et al. 2009; Modjaz et al. 2008b; Mazzali et al. 2008), which is shown by the thin arrow in the left panel (see Appendix 
A). 



X-ray transient or the subsequent tail (Soderberg et al. 
2008; Chevalier & Fransson 2008) while the later part 
(i > 4 days) is the SN component, powered by the decay 

of ^^Ni and ^^Co. 

The first part of the LC depends on the progenitor 
radius and radiation-hydrodynamics at outer layers, as 
well as Moj and i?K- To determine the global properties 
of the SN ejecta, we focus on the second, principal part, 
which depends on Mej, i?K and the amount of ejected 
^^Ni mass [M(^^Ni)]. The progenitor radius is discussed 
in §5.2. 

The LCs are calculated for the five explosion models 
presented in §2 (see Table 1). Our LTE, time-dependent 
radiative transfer code (Iwamoto et al. 2000) solves the 
Saha equation to obtain the ionization structure. Us- 
ing the calculated electron density, the Rossland mean 
opacity is calculated approximately by the empirical re- 
lation to the electron scattering opacity derived from the 
TOPS database (Magee et al. 1995, Deng et al. 2005). 
For the initial temperature structure of the SN ejecta, 
we use results of adiabatic hydrodynamic calculations at 
one day after the explosion. The hydrodynamics and the 
radiative transfer are not coupled. 

Asphericity of the ejecta of SN 2008D is suggested by 
the emission line profile in the spectrum at i = 109 days 
(Modjaz et al. 2008b). To include the possible effect of 
aspherical explosion, we modify the distribution of '"'^Ni 
from that derived from nucleosynthetic calculation. In 
hydrodynamic/nucleosynthetic calculations of aspherical 
explosion, more ^^Ni is mixed to the surface in the more 
aspherical cases (see e.g., Maeda et al. 2006, Tominaga 
2009). A constant mass fraction of ^^Ni is assumed below 
v^i, the outer boundary of ^^Ni distribution in velocity. 
The value of wni is determined so as to explain the rising 
part of the LC. The estimated i^ni is listed in Table 1. 
The resultant mass fraction of ^"^Ni is from 0.03 (HE4) 
to 0.01 (HE16). 

Figure 3 shows the calculated LCs compared with the 



observed LC. The model LCs of HE8, HEIO, and HE16 
reproduce the observed LC around the peak very well. 
The LCs of HE4 and HE6 tend to be narrower than the 
observations. At a later phase, the five LCs are all in 
good agreement with the observations. The steep decline 
in the calculated LCs at i < 4 days could be a relic of 
the shock-heated envelope, and radiation-hydrodynamics 
calculations are required to study this part. 

HE4 and HE6 need some enhancement of C in the 
He layer to reproduce the observed LC near the peak 
more nicely. The C-abundance in the He layer is poorly 
known because of the uncertainties involved in the C- 
production by convective 3 a-reaction in progenitor mod- 
els and those in the Rayleigh- Taylor instability at the 
He/C+0 interface during explosions, which tends to be 
stronger for lower mass He stars (Hachisu ct al. 1991). In 
view of these uncertainties, we include HE4 and HE6 in 
the further spectral analysis, rather than excluding them 
from the possible models. 

The timescale around the peak depends on both Afcj 

and i?K as oc k^^'^mJ E^ , where k is the optical 
opacity (Arnett 1982). Thus, for each model, a ki- 
netic energy can be specified so as to reproduce the ob- 
served timescale. The derived set of ejecta parameters 
are {M,;/Mq, Ek/W''^ erg) = (2.7, 1.1), (4.4, 3.7), (6.2, 
8.4), (7.7, 13.0) and (12.4, 26.5) for the case of HE4, 
HE6, HE8, HEIO and HE16, respectively. The ejected 
^^Ni mass is ^ O.O7M0 in all models. 

The model with M^j = 3 - 5Mq and Sk = 2 - 4 x lO^^ 
erg suggested by Soderberg et al. (2008) is close to our 
HE6 model, while the model of Mazzah et al. (2008), 
with Mej = 7Mq and Sk = 6 x 10^^ erg is close to 
our HE8 model. Model HE4 has the canonical explosion 
energy of core-collapse SNe (i.e., ~ 10^"'^ erg) while HEIO 
and HE 16 have the explosion energy of hypernovae (> 
10^2 erg). 

In all five models, the late evolution of the LC {t > 
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Fig. 4.— Spectrum of SN 2008D at t = 4.6 days from the X-ray 
transient (black line, Mazzali et al. 2008) compared with synthetic 
spectra (color lines). The spectra are shifted by 6.0, 4.5, 3.0, 1.5, 
0.0 from top to bottom. The model spectra are reddened with 
E{B — V) = 0.65 mag. From top to bottom, the synthetic spectra 
calculated with HE4, HE6, HE8, HEIO, and HE16 are shown. The 
red, blue and green lines show the synthetic spectra with oxygen 
mass fraction X{0) = 0.01, 0.1, and 0.5, respectively. Since the 
synthetic spectra with X{0) = 0.1 for more massive models than 
HE4 already show too strong O I line, the spectra with X{0) = 
0.5 are not shown for these models. 



200 days) is not very different, with a decline rate of 
^ 0.015 mag day~^ (right panel of Fig. 3). This decline 
is faster than the ^^Co decay rate (0.01 mag day~^, thin 
black line in Fig. 3) because some 7-rays escape without 
depositing energy in the SN ejecta at such late epochs. 
These models predict that the optical magnitude of SN 
2008D is ^ —10.5 mag (observed magnitude ~ 23.8 with 
no bolomctric correction) in 2008 October, i.e., ~ 300 
days after the explosion, when the SN can be observed 
again, and '^ —9.6 mag (~ 24.7 mag) at one year after 
the explosion (if dust does not form in the ejecta). 

4. OPTICAL SPECTRA 

In this section, the five models are tested against the 
observed spectra. Optical spectra have been shown by 
Soderberg et al. (2008), Malesani et al. (2009), Modjaz 
et al. (2008b) and Mazzah et al. (2008). We use the 
data set presented by Mazzali et al. (2008). The spectral 
sequence can be divided into three parts. At the earliest 
epochs (i < 4 days), the spectra are almost featureless 
^^. This is probably the result of shock heating (the first 
part of the LC). At 4 < t < 10 days, the spectra show 

^^ Two absorption features are identified around 4000 A in the 
spectra at t ~ 2 days (t = 1.77 days, Malesani et al. 2009; t = 1.84 
days, Modjaz et al. 2008b), while they are not seen in the spectra 
at t = 1.54 and 2.49 days presented by Mazzali et al. (2008). These 
absorptions might be due to more highly-ionized ions, such as C III, 



broad-line features. Around and after maximum (t > 10 
days), the spectrum shows strong He features as in Type 
lb SNe. The velocity of the He hues is - 9000 - 10, 000 
km s~^ (§4.3). We present spectral modelling at the SN 
dominated phase, i.e., i > 4 days. 

For spectral modelling, we use the one-dimensional 
Monte Carlo spectrum synthesis code (Mazzali & Lucy 
1993). The code assumes a spherically symmetric, 
sharply defined photosphere. Electron and line scatter- 
ing are taken into account. For line scattering, the ef- 
fect of line branching is included (Lucy 1999; Mazzali 
2000). The ionization structure is calculated with modi- 
fied nebular approximation as in Mazzali & Lucy (1993, 
see also Abbott & Lucy 1985). Although it is known that 
the non-thermal excitation is important for the He lines 
(Lucy 1991), non-thermal processes are not included in 
our analysis. Thus, we do not aim to obtain a good fit 
of the He lines. 

To determine the temperature structure, many photon 
packets are first traced above the photosphere with an 
assumed temperature structure. The Monte Carlo ray 
tracing gives the flux at each mesh and the temperature 
structure is then updated using the flux. This procedure 
is repeated until the temperature converges. Finally a 
model spectrum is obtained using a formal integral (Lucy 
1999). 

The input parameters of the code arc emergent lu- 
minosity (L), the position of the photosphere in veloc- 
ity (photospheric velocity, Uph), and element abundances 
(mass fractions) above the photosphere (i.e., in the SN 
atmosphere). Note that L and Wph do not depend much 
on the model parameters such as Mej and £'k. They are 
constrained by the absolute flux of the spectrum and the 
line velocities, respectively (and also by the relation of 
L ex I'pht^Jo'jj, where Tcff is the effective temperature of 
the spectrum). 

With the estimated luminosity and photospheric veloc- 
ity, mass fractions of elements are optimized. For sim- 
plicity, homogeneous abundances are assumed above the 
photosphere without using the results of nucleosynthctic 
calculations. We compare the derived abundances with 
those by nucleosynthetic calculations for the progenitor 
models. The goodness of the fit is judged by eyes be- 
cause of the complex dependences of the parameters and 
the difhculty in obtaining the perfect fit of the overall 
spectrum. 

4.1. Broad-Line Spectrum: At t = 4.6 Days 

We first perform model calculations for the spectrum 
at t = 4.6 days (Fig. 4). The spectrum shows broad-line 
features. 

4.1.1. Intermediate Mass Model HE8 

We use model HE8, the middle of our model sequence, 
as a fiducial case. A good agreement with the observed 
spectrum is obtained with tiph = 18, 500 km s^^ and log 
L (erg s""'^) = 41.7. Since this velocity is higher than the 

N III. and O III (Modjaz et al. 2008b; Quimby et al. 2007). We have 
investigated those lines by the Monte Carlo spectrum synthesis 
code, but we don't find a large contribution of these ions because 
ionizations by the photospheric radiation only is not enough for 
the strong contributions of such ions, as noted by Modjaz et al. 
(2008b). 
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He line velocities observed later phases {^ 9000-10,000 
km s^^), the photosphere at this epoch is expected to be 
located in the He-rich layer. 

Figure 4 shows a comparison of the observed and syn- 
thetic spectrum. The spectrum has P-Cygni profiles of 
O I, Na I, Ca ii, Ti ii, Cr ii, and Fe ii lines. The line 
at 6000 A is identified as Si II. The contribution of the 
high velocity Ha is quite small, which is discussed in 
Appendix B. The spectrum at wavelengths bluer than 
5500A is dominated by Ti ii, Cr ii, and Fe ii lines. Given 
the uncertainty in the metal abundances in outer layers, 
reflecting the uncertainty of the explosion mechanism or 
the degree of mixing, these features can be fitted using 
the optimal value of the metal abundances. 

In contrast to the heavy, synthesized elements, the oxy- 
gen abundance cannot be totally parameterized because 
the majority of oxygen is synthesized during the evolu- 
tion of the progenitor star. The red and blue line shows 
synthetic spectra with oxygen abundance X(O)=0.01 
and 0.1, respectively. In these models, the abundance 
of He is X(Hc)~0.8 and 0.7, respectively. The spectrum 
with X(O)=0.01 (red) gives a good match with the ob- 
served O iA7774 line around 7400A, while the O i line in 
the model spectrum with X(O)=0.1 (blue) is too strong. 
Since the oxygen abundance in the He-rich layer is of the 
order of 10~^ almost irrespective of evolutionary models, 
this is consistent with the fact that the photosphere is 
located in the He-rich layer. 

In the observed spectrum, the O i and Ca ii IR triplet 
are blended at 7000 - 8500 A while they are separated in 
the synthetic spectra. This is caused by the insufficient 
Ca II absorption in the model at the very high velocity 
layer with v ~ 0.1c. The ejecta mass at v > 0.1c in HE8 
is O.O4M0 (Table 1), which is consistent with that in the 
model presented by Mazzali et al. (2008). Note that the 
mass at w > 0.1c is much smaller than in the model for 
SN 1998bw (- 1.5Mq, CO138E50 in Nakamura et al. 
2001a, see also Fig. 2). 

4.1.2. Massive Models HEIO and HE16 

Next, we use the more massive models. For model 
HEIO, the O I line in the model with X(O)=0.01 (red) 
seems to be strong, but we can obtain a good fit with 
slightly smaller oxygen abundance. For model HE16, the 
strength of the O i line with X(O)=0.01 (red) is similar 
to that in HEIO. In these massive models, the O I feature 
is too strong in the model spectra with X(O)=0.1 (blue). 
This is consistent with the fact that the photosphere (v = 
18, 500 km s^^) is located in the He-rich layer. 

4.1.3. Less Massive Models HE6 and HE4 

Finally, we use less massive models. For the less mas- 
sive model HE6, X(O)=0.01 gives a reasonable fit to the 
O I line. In contrast, X(O)=0.1 yields too strong a line. 
This is a similar behavior to the more massive models, 
and implies that the photosphere is located in the He-rich 
layer. 

For HE4, the synthetic spectra with X(O)=0.01 and 
0.1 do not give a strong enough O i absorption (red 
and blue lines). To explain the observed absorption, 
X(O)=0.5 is required (green line) because of the low 
density at the outer layer of HE4 (Fig. 2). This requires 
that the layer at u = 18,500 km s^^ should already be 
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Fig. 5.— Spectrum of SN 2008D at t = 32.4 days from the X-ray 
transient (black line, Mazzali et al. 2008) compared with synthetic 
spectra (color lines). The spectra are shifted by 20.0, 15.0, 10.0, 
5.0, 0.0 from top to bottom. The model spectra are reddened with 
E{B — V) = 0.65 mag. From top to bottom, the synthetic spectra 
calculated with HE4, HE6, HE8, HEIO, and HE16 are shown. The 
red and blue lines show the synthetic spectra with an oxygen mass 
fraction of X{0) = 0.01 and 0.8, respectively. Dashed, vertical 
lines show the position of the He lines (He I 5876, 6678, 7065) 
blueshifted with v = 10,000 km s~^. Note that the He lines are 
not treated in the code. 



0-rich, which is clearly inconsistent with the observed 
He line velocity {v ~ 9000-10,000 km s^^). Therefore, 
HE4 is not likely to be a viable model for SN 2008D. 

4.2. Type lb Spectrum: At t = 32.4 Days 

At < = 32.4 days, the observed spectrum shows typi- 
cal Type lb features. The overall features are fitted well 
with Wph = 7500km s"^ and log L (erg s'^) = 42.1 (Fig. 
5). The observed Fe lines at 4500 - 5000 A arc too nar- 
row to be reproduced by the massive models (HE8, HEIO 
and HE16). This is caused by our crude assumption of 
a homogeneous abundance distribution, which is not ap- 
propriate for synthesized elements such as Fe in the inner 
layer. The spectra might be improved using a stratified 
abundance distribution or non-spherical models (Tanaka 
et al. 2007). 

4.2.1. Intermediate Mass Model HE8 

The model spectrum with X(O)=0.01, as assumed for 
the spectrum at t=4.6 days, is shown in red line. The 
synthetic O i line at 7500 A is slightly weaker than 
the observation. A value Ar(O)=0.8 yields a reasonably 
strong O I line (blue) , which implies that the photosphere 
{v — 7500 km s^^) is not located in the He-rich layer, 
consistent with vn^ = 10500 km s~^ (velocity at the bot- 
tom of the He layer) of HE8. 
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Fig. 6.— Spectral evolution of SN 2008D (black, Mazzali et 
al. 2008) compared with the sequence of the synthetic spectra 
computed with HE8 (red). The epoch in the figure shows the days 
from the X-ray transient. The values in the parenthesis shows the 
mass fraction of oxygen assumed in the calculation. The spectra 
are shifted by 16.0, 14.0, 12.0, 10.0, 8.0, 6.0, 4.0, 2.0, 0.0 from top 
to bottom. The model spectra are reddened with E{B — V) = 0.65 
mag. Dashed, vertical lines show the position of the He lines (He I 
5876, 6678, 7065) blueshiftcd with v = 10,000 km s'^. Note that 
the He lines are not treated in the code. 



4.2.2. Massive Models HEIO and HE16 

The synthetic spectra calculated using HEIO and HE16 
show the stronger O i line than HE8. The spectrum with 
X(O)=0.01 (red) gives a slightly weaker O i line than 
in the observation, while the spectrum with X(O)=:0.8 
(blue) yields a sufficiently strong line. Although the syn- 
thetic O I line with X(O)=0.8 is too strong especially at 
high velocity (i.e., at bluer wavelength), this is caused by 
the assumption of the homogeneous abundance distribu- 
tion. Thus, near the photosphere, a high mass fraction 
of O is preferred. This is consistent with the high whc 
of these models (12, 500 and 17, 500 for HEIO and HE16, 
respectively) . 

However, the observed He line velocities {v ^ 9000- 
10,000 km s~i) suggest that the layer at w ~ 10,000 
km s~^ is still He-rich. This is inconsistent with the 
high uhc in HEIO and HE16, requiring that the layers at 
V ~ 10,000 km s-i be 0-rich. 

4.2.3. Less Massive Models HE6 and HE4 

The synthetic spectra using HE6 also have similar 
trend with those of HE8. The spectrum with X(O)=0.8 
(blue) gives a reasonable fit to the O i line. Although 
the low velocity at the bottom of the He layer in HE6 
(uhc = 6700 km s~^) suggests that the photosphere at 
this epoch {v = 7500 km s~^) is still in the He layer, this 
small difference is within the uncertainty of whc caused 
by the variation of the He layer mass depending on evo- 
lutionary models. 
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Fig. 7. — Time evolution of photospheric velocity calculated 
by the LC code (color lines) and that derived from the spectral 
modelling (filled circles). Open black symbols show the line velocity 
of He lines. The horizontal lines show i>He for each model (Table 
1). Note that ?;ho can be varied depending on evolutionary models 
with different masses of the He layer. 



For HE4, the O i absorption is reproduced with X{0)= 
0.8. However, the very low fHc of HE4 (< 3500 km s~^) is 
not consistent with the fact that the spectrum model re- 
quires the 0-dominated photosphere at w = 7500 km s~^. 

4.3. Velocity Evolution 

Using the fiducial model HE8, we calculate the spec- 
tral evolution (Fig. 6). The values in the parenthesis 
is the oxygen mass fraction adopted in the fitting. We 
find that a higher oxygen mass fraction is preferred for 
the later spectra. Although a homogeneous mass frac- 
tion is assumed in the calculation, the photospheric posi- 
tion seems to transit from the He-rich layer to the 0-rich 
layer around t = 25.5 days. Thus, the boundary between 
the He-rich and 0-rich layers is located near v = 7800 
km s~^. 

Figure 7 shows the photospheric velocities derived from 
the spectral modelling (filled black circles), which does 
not depend much on the model parameters. In Figure 7, 
the photospheric velocities obtained from the synthetic 
LCs (§3) are also shown (solid lines). The photospheric 
velocities for HE6 and HE8 are close to the values de- 
rived from spectral modelling. However, it should be 
noted that the photospheric velocities obtained from the 
LC models are only approximate because the LC model 
assumes LTE and does not fully take into account the 
contribution of the line opacity. Thus, uncertainty of a 
few thousand km s~^ is expected. Nevertheless, the de- 
creasing trend of the photospheric velocity derived from 
the spectral modelling is reproduced by our LC calcula- 
tions because (1) we use the hydrodynamic models (Fig. 
2) having a decreasing density structure toward the outer 
layers and (2) we solve the ionization in the ejecta. and 
thus, the opacity is time-dependent. 

In Figure 7, the Doppler velocities of three He i lines 
measured at the absorption minimum (open symbols) are 
also shown. Malesani et al. (2009) and Modjaz et al. 
(2008b) show the subsequent spectral evolution, and the 
He line velocity declines slowly to u ~ 9000 km s~^. 
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Fig. 8. — Kinetic energy of the explosion (upper) and the ejected 
^^Ni mass as a function of the estimated main-sequence mass of 
the progenitors for several core-collapse SNe. The parameters are 
listed in Table 2 with references. The progenitor mass of SNe 
shown in the figure is estimated based on the Mms-Mc relation by 
Sugimoto & Nomoto (1980; used in Nomoto & Hashimoto 1988) 
as in this paper. 



The horizontal Hnes in Figure 7 mark the velocity at 
the bottom of the He layer for the five models (whc, see 
Table 1). In HE 10 and HE16, whc is too high compared 
with the observed velocities (§4.2.2). Also in HE8, it 
may be higher than the minimum of the observed He line 
velocity {v "^ 9000 km s~^). The lower vue in HE4 and 
HE6 cannot be excluded from the observed line velocities. 
But the spectral modelling shows that the layer at w ~ 
7500 km s^^ is not He-rich (§4.2). This is inconsistent 
with the very low vjje in HE4. It must be cautioned that 
uhc is affected by the mass of the He layer (i.e., by the 
choice of evolutionary models) . 

5. DISCUSSION 

5.1. Optimal Model for SN 2008D 

For the five He star progenitor models, we calculate hy- 
drodynamics of the explosions and explosive nucleosyn- 
thesis. To reproduce the observed LC, we obtain the 
possible set of the mass and kinetic energy of the ejecta: 
(Mej/Mo, Sk/IO^i erg) = (2.7, 1.1), (4.4, 3.7), (6.2, 8.4), 
(7.7, 13.0) and (12.4, 26.5) for HE4, HE6, HE8, HEIO 
and HE16, respectively. These five models are tested 



against the optical spectra. 

Model HE4 has many difficulties in reproducing the 
observed spectra. At early epochs, the calculated O i 
line is too weak because of the too small oxygen mass 
in the He-rich layer. At later epoch, the model spec- 
trum suggests that the photospheric layer at w = 7500 
km s~^ is 0-rich, which is not consistent with the ex- 
plosion model that has He-rich or He-0 mixed layers at 
V > 3000kms-i. 

Model HE6 can reproduce the observed spectra well. 
The evolution of the photospheric velocity calculated 
with HE6 is in reasonable agreement with the velocities 
derived from the spectral modelling (Fig. 7). The spec- 
tral model at t = 32.4 suggests that the layer sA v ^ 7500 
km s~^ is not the He-rich layer, while the slightly lower 
v-ne of HE6 (6700 km s^^) implies that the photosphere 
at this epoch is He-rich. 

Model HE8 is reasonably consistent with all the as- 
pects studied in this paper At all epochs, the optical 
spectra can be explained with a reasonable abundance 
distribution, and the calculated photospheric velocities 
are consistent with those derived from spectrum synthe- 
sis. However, the velocity at the bottom of the He layer 
(t^He) in HE8 is slightly higher than the observed He line 
velocities. 

Model HEIO and HE16 reproduce the early and later 
spectra reasonably well. However, these models predict 
too high photospheric velocity (Fig. 7). In addition, 
the velocities at the bottom the He layer (who = 12, 500 
and 17, 500 km s^^for HEIO and HE16, respectively), are 
not consistent with the observed line velocity {v ~ 9000 
- 10,000kms-i). 

In summary HE4, HEIO and HE 16 are not consistent 
with SN 2008D. Both HE6 and HE8 have a smaU incon- 
sistency related to the boundary between the He-rich and 
0-rich layers. It seems that a model between HE6 and 
HE8 may be preferable. However, since there is uncer- 
tainty in WHe in our model set, depending on the mass of 
the He layers, we include both HE6 and HE8 as possible 
models. 

We conclude that the progenitor star of SN 2008D has 
a He core mass M^ = 6 — 8M0 prior to the explosion. 
This corresponds to a main-sequence mass of Mms = 
20 — 25Mq under the Myis-Ma relation by Sugimoto & 
Nomoto (1980; used in Nomoto & Hashimoto 1988). We 
find that SN 2008D is an explosion with Moj = 5.3 ± 
l.OMg and E]i = 6.0 ± 2.5 x lO^^ erg. The mass of 
the central remnant is 1.6 — l.SM©, which is near the 
boundary mass between the neutron star and the black 
hole. Note that the error bars only refiect the uncertainty 
of the LC and spectral modelling. Possible additional 
uncertainties of the parameters are discussed below. 

Distance and reddening: since the distance to the 
host galaxy and the reddening toward the SN include 
some uncertainties, the ejected ^^Ni mass could also con- 
tain ^ 20% uncertainties. However, Moj and Mcut are 
not affected because these values are much larger than 
the ejected ^^Ni mass. Thus, the estimated core mass 
and progenitor mass are not largely affected by the un- 
certainty of the distance and the reddening. 

Asphericity of the explosion: possible effects on 
the estimate of Aioj and -Bk from asphericity of the ejecta 
are of interest. These effects were studied for SN 1998bw 
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TABLE 2 
Parameters of Supernovae 



SN (Type) 


Mci ^ 


Ek"^ 


M(56Ni) c 


Mms ^ 


Rcfs. 


SN 1987A (11 pec) 


14.7 


1.1 ±0.3 


0.07 


20 ±2 


1, 2 


SN 1993J (lib) 


3.2 ±0.3 


1.1±0.1 


0.08 ±0.03 


13.5 ±1.5 


3 


SN 19941 (Ic) 


1.05 ±0.15 


1.0 ±0.2 


0.07 ±0.03 


14 ±1 


4, 5 


SN 1997D (II) 


~24 


~0.4 


0.0045 ± 0.0035 


30 ±10 


6 


SN 1997ef (Ic) 


8.6 ±1 


12.75 ± 4.75 


0.15 ±0.03 


32.5 ±2.5 


7, 8 


SN 1998bw (Ic) 


10.4 ±1 


40 ± 10 ° 


0.43 ±0.05 


40 ±5 


9, 10 


SN 1999br (II) 


~ 14 


~ 0.6 


0.002 


25 ±9 


11 


SN 2002ap (Ic) 


3.25 ±0.75 


4±1 


0.075 ± 0.005 


22.5 ±2.5 


12 


SN 2003dh (Ic) 


7±3 


35 ±15 


0.4 ±0.125 


32.5 ±7.5 


13, 14 


SN 20031w (Ic) 


~13 


55 ±5 


0.55 ±0.05 


45 ±5 


15 


SN 2005bf (lb pec) 


6.5 ±0.5 


1.25 ±0.25 


0.04 ±0.02 f 


25 ±2 


16, 17 


SN 2006aj (Ic) 


1.8 ±0.8 


2.0 ±1.0 


0.21 ±0.04 


20 ±2 


18 


SN 2008D (lb) 


5.3 ±1.0 


6.0 ±2.5 


0.07 ±0.005 


22.5 ±2.5 


this work 



References. — (1) Shigoyama & Nomoto (1990), (2) Blinnikov et al. (2000), (3) Shigeyama et al. (1994), (4) Iwamoto et al. (1994), 
(5) Sauer et al. (2006), (6) Turatto ct al. (1998), (7) Iwamoto et al. (2000), (8) Mazzali, Iwamoto & Nomoto (2000), (9) Nakamura et al. 
(2001a), (10) Iwamoto et al. (1998), (11) Zampieri et al. (2003), (12) Mazzali et al. (2002), (13) Mazzali et al. (2003), (14) Deng et al. 
(2005), (15) Mazzali et al. (2006a), (16) Tominaga ct al. (2005), (17) Maeda ct al. (2007), (18) Mazzali et al. (2006b) 

''The mass of the SN ejecta (Mq) 

''The kinetic energy of the SN ejecta (10^^ erg) 

<^The mass of ejected 56Ni (Mq) 

"* Estimated main-sequence mass (Mq) 

"^i^K = 20 X 10^^ erg is derived from the modelling with a multi-dimensional model (in the polar-viewed case, Maeda et al. 2006; Tanaka 
et al. 2007). 

^The mass of ^^Ni is derived from the late time observation (Maeda et al. 2007). The early observations suggest M(^''Ni)= O.3M0 
(Tominaga ct al. 2005; Folatelli et al. 2006). 



associated with GRB 980425 by Maeda et al. (2006) and 
Tanaka et al. (2007). They found that the kinetic energy 
can be smaller by a factor of < 2 in the on-axis case of 
highly aspherical explosion than in the spherical model. 
There is little effect in the off-axis case. 

Modjaz et al. (2008b) presented the spectrum at 
t = 109 days and suggested the asphericity of SN 2008D 
by the doubly-peaked emission profile of the O i line. 
Such a profile of the O i line has been interpreted as 
an off-axis line-of-sight in the axisymmetric explosion 
(Maeda et al. 2002; MazzaH et al. 2005; Maeda et al. 
2008; Modjaz et al. 2008a). Thus, the estimate by the 
modelling under spherical symmetry may not be largely 
changed even for the aspherical models. The quantitative 
discussion should wait for the later spectra, the detailed 
modelling of the line profile, and determination of the 
degree of asphericity and the line-of-sight. 

The effects of asphericity, especially aspherical mass 
ejection and fallback, are also important to determine the 
relation between the ejected ^^Ni mass and the remnant 
mass. The remnant mass in this paper is determined 
to eject the optimal amount of ^^Ni by one-dimensional 
hydrodynamic/nucleosynthetic calculations. However, 
since the remnant mass could be either larger or smaller 
depending on the asphericity and details of the explosion 
mechanism, the estimate by one-dimensional calculations 
is a reasonable approximation. 

Possible presence of hydrogen: Soderberg et al. 
(2008) identified the high velocity Ha for the absorption 
line at 6150 A. If the mass of the H layer is not negligible, 
it might affect the core mass, which we estimate by as- 
suming non-existence of H, i.e., a bare He core. However, 
we find a large mass of the H layer is inconsistent with 
the spectrum at i = 4.6 days (Appendix B). The mass of 
the H layer is smaller than 5 x 10~^Mq, and thus, there 



is no effect on the parameters. 

Evolutionary models: the estimate of the main- 
sequence mass uses the approximate Mms-Mqi relation 
by Sugimoto & Nomoto (1980; Eq. 4.1), which is used 
in Nomoto & Hashimoto (1988). The Mms-Mc, rela- 
tion of several evolutionary models are shown in Fig- 
ure 1. The systematic differences in this relation for 
^MS < 3OM0 (Ma < IOMq) may stem from the differ- 
ences in the treatment of convection, mass loss, rotation, 
and binary effects (§2.1). Thus, we should keep in mind 
that the main-sequence mass is subject to systematic un- 
certainties of 3 — 5Mq (Fig. 1). Note that our estimate 
of the He core mass depends only on the estimates of 
Mej and -Ek from hydrodynamic/nucleosynthetic calcu- 
lations, and thus, our determination of the He core mass 
is not affected by the variety of the evolutionary models. 

5.2. Comparison with Previous Works 

Soderberg et al. (2008) have estimated the parameters 
of the ejecta as Mej = 3 - 5Mq and Ek = 2 - 4 x 10^^ 
ergs, which are smaller than those derived in this pa- 
per. The difference seems to stem from their assumptions 
of the homogeneous sphere and time-independent opac- 
ity. These assumptions lead an almost time-independent 
photospheric velocity, which is not the case in SNe. 
Especially for SN 2008D, the very early spectra show 
the broad-line features, and the photospheric velocity at 
t ^ 5 days after the X-ray transient is almost twice as 
high as the velocity around maximum. 

Adopting the cooling envelope model by Waxman et 
al. (2007) to the blackbody temperature and the radius 
at < < 4 days, Soderberg et al. (2008) estimated a pro- 
genitor radius to be R^ ~ li?© with E{B — V) = 0.61 
mag, Moj = 5Mq and Ek ^ 2 x 10^^ erg. Modjaz et al. 
(2008b) also derived a similar value, i?* = 1.1 ± O.46i?0 
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Fig. 9. — Spectral comparison among SNe 1993J (lib; Barbon 
et al. 1995), 1999ex (Ib/c; Hamuy et al. 2002), 2008D (lb) and 
2005bf (peculiar lb; Anupama et al. 2005). The epoch for SNc 
1993J, 1999ex, and 2005bf is given in the estimated days from 
the explosion. Dashed, vertical lines show the position of the He 
lines (He I 5876, 6678, 7065) blueshifted with v = 10,000 km s'^ 
The spectra except for SN 2008D are taken from the SUSPECT 
database. 



with E{B — V) — 0.6 mag and the same Moj and Ey^ 
with Soderberg et al. (2008). If Mgj and Ey^ derived in 
this paper are adopted, the estimated radius is ~ 80% 
of their estimate. This is marginally consistent with the 
radius of model HE8 while it is smaller than HE6. In 
this sense, model IIE8 seems to be more self-consistent. 
It must be noted, however, that Chevalier & Fransson 
(2008) derived a larger radius, i?* ~ 9i?0 by using the 
model by Chevalier (1992, and using the blackbody tem- 
perature and the radius presented by Soderberg et al. 
2008). 

Mazzali et al. (2008) estimated the ejecta parame- 
ters by modelling the bolometric LC and optical spec- 
tra. Their largest assumption is that a central rem- 
nant as massive as ~ SA/q is implicitly assumed, which 
leads a massive He core mass (~ IOMq), and thus, 
a massive progenitor mass (~ 30Mq). Our hydrody- 
namic/nucleosynthetic calculations show that a smaller 
central remnant is preferred (^ 1.6 — 1.8Mq, Table 1). 

5.3. SN 2008D in the Context of Type Ib/c Supernovae 

Figure 8 shows the kinetic energy of the ejecta and the 
ejected ^^Ni mass as a function of the estimated main- 
sequence mass for several core-collapse SNe (see, e.g., 
Nomoto et al. 2007). The parameters shown in Figure 
8 are also listed in Table 2. SN 2008D is shown by a 
red circle in Figure 8. The ejecta parameters for other 
SNe shown in Figure 8 and Table 2 are derived from 
one-dimensional modelling as in this paper. Although 
there is a systematic uncertainty in the progenitor mass 
(Fig. 1), the progenitor mass of SNe shown in Figure 8 
is estimated based on the Mms-M^ relation by Sugimoto 



& Nomoto (1980; used in Nomoto & Hashimoto 1988) as 
in this paper. Thus, the relative position of SNe in the 
plots are robust. 

The main-sequence mass of the progenitor of SN 2008D 
is estimated to be between normal SNe and GRB-SNe 
(or hypcrnovae). The kinetic energy of SN 2008D is also 
intermediate. Thus, SN 2008D is located between the 
normal SNe and the "hypernovae branch" in the -Ek — 
Mms diagram (upper panel of Fig. 8). The ejected ^^Ni 
mass in SN 2008D (- 0.07Mq) is similar to the ^^m 
masses ejected by normal SNe and much smaller than 
those in GRB-SNe. 

Figure 9 compares the spectra of SNe 1993J (lib, Bar- 
bon et al. 1995), 1999ex (Ib/c, Hamuy et al. 2002), 
2005bf (lb, Anupama et al. 2005, Tominaga ct al. 2005, 
Folatclh ct al. 2006), and 2008D (lb). The epoch for SNe 
1993 J, 1999ex, and 2005bf is given in the estimated days 
from the explosion. The explosion epoch is uncertain up 
to - 15 days in SN 2005bf (FolateUi et al. 2006) while it 
is well constrained in SNe 1993J and 1999ex (< 2 days, 
Wheeler et al. 1993; Hamuy et al. 2002). 

The spectra of SN 2008D and SN 1999ex are very sim- 
ilar (Valcnti et al. 2008b), while SN 2005bf has a lower 
He velocities. Although the epoch of SN 2005bf is uncer- 
tain, the He line velocities in SN 2005bf is always lower 
than 8000 km s'^ (Tominaga et al. 2005). The He hues 
in SN 1993J are very weak at this epoch. The Fe features 
at 4500-5000A are similar in these four SNe, but those 
in SN 2005bf are narrower. 

Malesani et al. (2009) suggested that the bolometric 
LCs of SNc 1999ex and 2008D are similar. If it is the 
case (although some discrepancy is shown by Modjaz et 
al. 2008b). the similarity in both the LC and the spectra 
suggests that SN 1999cx is located close to SN 2008D in 
the Ek - Mms and MC^^Ni) - Mus diagrams. 

Comparison with other Type lb SNe shown in Figure 
8 is possible only for SN 2005bf although SN 2005bf is 
a very peculiar SN that shows a double peak LC with 
a very steep decline after the maximum, and increasing 
He line velocities (Anupama et al. 2005; Tominaga et al. 
2005; FolateUi et al. 2006; Maeda et al. 2007). The LC 
of SN 2005bf is broader than that of SN 2008D, while 
the expansion velocity of SN 2005bf is lower than that 
of SN 2008D. These facts suggest that SN 2005bf is the 
explosion with lower Ek/Mcj ratio (Table 2). 

Malesani et al. (2009) also pointed the similarity of 
the LCs of SNe 1993J and 2008D. But the expansion 
velocity is higher in SN 2008D (see, e.g., Barbon et al. 
1995; Prabhu et al. 1995). Thus, both the mass and the 
kinetic energy of the ejecta are expected to be smaller in 
SN 1993J. In fact, SN 1993J is explained by the explosion 
of a 4:Mq He core with a small mass H-rich envelope 
(Nomoto et al. 1993; Shigeyama et al. 1994; Woosley et 
al. 1994). 

6. CONCLUSIONS 

We presented a theoretical model for SN 2008D asso- 
ciated with the luminous X-ray transient 080109. Based 
on the progenitor models, hydrodynamics and explosive 
nucleosynthesis are calculated. Using the explosion mod- 
els, radiative transfer calculations are performed. These 
models arc tested against the bolometric LC and optical 
spectra. This is the first detailed model calculation for 
Type lb SN that is discovered shortly after the explosion. 



Type lb Supernova 2008D 



11 



We find that SN 2008D is a more energetic explosion 
than normal core-collapse SNe. We estimate that the 
ejecta mass is Mej = 5.3 ± I.OMq and the total kinetic 
energy of Ek = 6.0 ± 2.5 x lO^i erg. The ejected ^^Ni 
mass is ~ O.OTMq. To eject the optimal amount of ^^Ni, 
the mass of the central remnant is estimated to be 1.6 — 
I.SMq. The error bars include only the uncertainty of 
the LC and spectral modelling. 

Summing up the above masses, it is concluded that 
the progenitor star of SN 2008D has a 6 — 8Mq He core 
prior to the explosion. There is essentially no H enve- 
lope with the upper limit of 5 x 10~'^Mq. Thus, the 
corresponding main-sequence mass of the progenitor is 
Mms = 20 — 25M0 under the AIus-Ma relation by Sug- 
imoto & Nomoto (1980, used in Nomoto & Hashimoto 
1988). We note that there exist additional systematic 
uncertainties in this relation due to convection, mass 
loss, rotation, and binary effects. Our estimates of these 
masses and energy suggest that SN 2008D is near the bor- 



der between neutron star-forming and black hole-forming 
SNe, and has properties intermediate between those of 
normal SNe and hypernovae associated with gamma-ray 
bursts. 
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Fig. B1. — Left: Optical spectrum at t = 18.5 days (Mazzali ct al. 2008) compared with the synthetic spectra. The rod and blue lines 
show the synthetic spectra with the models where the Si mass fraction is twice as large as the solar abundance (best fit) and the same as 
the solar abundance, respectively. The element abundances are assumed to be homogeneous in these models. The green and magenta lines 
show the synthetic spectra with X(H)=0.8 and 0.4 at ?; > 18000 km s~^, respectively. The corresponding mass of H is 0.4 and 0.2 Mq, 
respectively. The homogeneous, solar abundance of Si is assumed in these models. Right: The optical spectrum at t = 4.6 days (Mazzali 
et al. 2008) compared with the synthetic spectra. The red line shows the model with the solar abundance of Si, which gives a good fit to 
the observed spectrum. The green, blue, magenta and cyan lines show the models with X(H) = 0.4,0.05,0.005 and 0.001 at v > 18000 
km s~^, respectively. The corresponding mass of H is 0.4, 0.025, 0.0025 and 0.0005 Mq, respectively. 



APPENDIX 
A. CONSTRUCTION OF BOLOMETRIC LIGHT CURVE 

The bolomctric LC shown in this paper was eonstructed by using optical data taken by the MAGNUM teleseope 
(Yoshii 2002; Yoshii, Kobayashi & Minezaki 2003), the Himalayan Chandra Telescope, and Swift UVOT (U-band, 
Soderberg et al. 2008), and also NIR data taken by the MAGNUM telescope. 

The bolometric luminosity was derived integrating the flux from U (with the edge of 9.68 x lO^** Hz) to K (1.00 x 10^** 
Hz) band. The photometric points are interpolated by the third order natural spline. If the data point of a certain 
band is not available, we use linear interpolation of the magnitude. 

The derived bolometric LC can be compared with that by Soderberg et al. (2008), Malesani et al. (2009), Modjaz 
et al. (2008b, U-Ks integration) and Mazzali et al. (2008). Although the scatter up to 0.4 mag is found among the 
LCs by direct comparison, it is caused mainly by the difference in the assumed distance and reddening. 

If the same distance and reddening are used (here we approximately correct the difference in the bolometric magnitude 
AMboi caused by the difference in the assumed reddening by AMboi = Rv^E{B — V), where Rv~3.1 and AE{B~V) 
is the difference in the assumed color excess), the LCs around/after the maximum is consistent among those by this 
paper, Soderberg et al. (2008), Modjaz et al. (2008b) and MazzaH et al. (2008) within 0.1 mag, while the LC by 
Malesani et al. (2009) is fainter by 0.2-0.5 mag. For the pre-maximum epochs, the LCs by these papers are consistent 
within 0.2 mag except that the magnitude at t = 4 days by this paper (shown by the arrow in the left panel of Fig. 3) 
is brighter than other ones by 0.25 mag. 

Since the scatter in the maximum luminosity among the papers is up to 0.2 mag, it causes the uncertainty of the 
ejected ^®Ni mass up to ^ 20%. However, this uncertainty does not affect our determination of the ejecta mass because 
the change in the ^^Ni mass (and mass cut) is negligible compared to the ejecta mass (§5). In addition, the time scale 
of the bolometric LC around the maximum is reasonably consistent among the papers, the kinetic energy of the ejecta 
is also not affected. 



B. NONEXISTENCE OF THE HYDROGEN LAYERS 

Soderberg et al. (2008) identified the high velocity (HV) Ha line for the absorption feature around 6150 A in the 
spectra around maximum. It is blended with the strong Si li line, and discrimination is not easy (e.g.. Branch et al. 
2006; Elmhamdi et al. 2006). The presence of H is important to specify the properties of the progenitor star just prior 
to the explosion. In addition, if the H layer is present, the estimate of A/cj and Ek may be affected since we have used 
bare He stars for the LC and spectral modelling. 

First, we test the presence of H in the spectrum around maximum using model HE8. The left panel of Figure Bl 
shows the comparison between the observed spectrum at t = 18.5 days (Mazzali et al. 2008) and synthetic spectra. 
The photospheric velocity at this epoch is 9000 km s^^ (Fig. 7). If the absorption at 6150 A is Si ii A6355, the 
Doppler velocity of the absorption at 6150 A is 9300 km s"^, which is well consistent with the photospheric velocity. 
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The red line shows the best fit model that includes Si twice as large as the solar abundance. The absorption is slightly 
shallower in the model with solar abundance Si (blue). Since the abundance twice as large as the solar abundance is 
reasonable for the middle layers of the ejecta, the HV H i is not necessarily required. 

However, this does not exclude the possibility of the presence of H at the outer layers. If the absorption at 6150 A 
is Ha, the Doppler velocity is 18500 km s^^. To test the presence of the H at such high velocity layers, we calculate 
model spectra by replacing He at w > 18000 km s~^ with H. The green and magenta lines show the models with 
X(H)=0.8 and 0.4 at z; > 18000 km s~^, respectively. The corresponding mass of H is 0.4 and 0.2 Mq, respectively. 
The models also include the solar abundance of Si at at u > 9000 km s~^. While the model with X(H)=0.8 gives a 
too strong absorption, the model with X(H)=0.4 agree with the observed spectrum. Thus, the presence of 0.2 Mq of 
H cannot be denied from the spectrum around the maximum. 

Next, we perform the similar tests using the very early spectrum. The right panel of Figure Bl shows the comparison 
between the observed spectrum at t = 4.6 days (Mazzali et al. 2008) and the synthetic spectra. The red line shows the 
best fit model, which include the solar abundance of Si. The blue line shows the model that have the solar abundance 
of Si and X(H)=0.4 at w > 18000 km s^^. Although this model gives a reasonable fit to the spectrum at i=18.5 
days (left panel of Fig. Bl), it shows too strong Ha and H/? at t = 4.6 days (the lack of H/? has been pointed out 
by Malessani et al. 2009). We get the stronger line at earlier epochs because the density at the high velocity layers 
{v = 18000 km s~^) becomes lower with time, and the line forming there is more effective at earlier epochs. 

The green, magenta and cyan lines show the models with smaller mass fraction of H, X(H)=0.05, 0.005, and 0.001, 
respectively. The corresponding mass of H is 0.025, 0.0025 and 0.0005 Mq. The models with X(H) = 0.05 and 0.005 
(green and magenta) still shows too strong H i lines. With X(H) — 0.001, the Ha line has little effect on the absorption 
at 6150 A although the model spectrum still has a sharp absorption of the HV Ha. 

If we use model HE6, the mass at the outer layers is smaller. Thus, we conclude that the mass of H is smaller than 
5 X 10-^ Mq. 



